Bending crystals: Emergence of fractal dislocation structures 
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We provide a minimal continuum model for mesoscale plasticity, explaining the cellular disloca- 
tion structures observed in deformed crystals. Our dislocation density tensor evolves from random, 
smooth initial conditions to form self-similar structures strikingly similar to those seen experimen- 
tally - reproducing both the fractal morphologies and some features of the scaling of cell sizes and 
misorientations analyzed experimentally. Our model provides a framework for understanding emer- 
gent dislocation structures on the mesoscale, a bridge across a computationally demanding mesoscale 
gap in the multiscale modeling program, and a new example of self-similar structure formation in 
non-equilibrium systems. 
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Structural engineering materials have a bewildering va- 
riety of microstructures, which are often controlled by 
deformation and annealing during the formation process. 
An imposed distortion generates a complex morphology 
even for a single crystal of a pure material - polycrys- 
tallinc grains form at high temperature where disloca- 
tion climb allows for polygonization, cell structures form 
at low temperatures when climb is forbidden. Cell walls 
(Fig. Hhjd) are distinct from grain boundaries in that 
they have smaller misorientations, different origins, are 
morphologically fuzzier, and the cells refine (get smaller) 
under shear. Experiments differ in characterizing the 
cell structures; some show convincing evidence of frac- 
tality [3-i] with structure on all length scales (Fig. [1}:), 
while others show structures with a single characteristic 
scale setting their cell size and cell wall misorientation 
distributions d-Q (Fig.[TU)- 

Dislocation avalanches 0], size-dependent hardness 
(smaller is stronger) [H and cellular structures 0, Q all 
emerge from collective dislocation interactions on the mi- 
cron scale. We expect that these mesoscale phenomena 
should be captured by an appropriate continuum the- 
ory of dislocation dynamics. Computationally, such a 
theory is crucial for multiscale modeling, as atomistic 
and discrete dislocation simulations are challenging on 
these scales of length and strain. Here we present a min- 
imal model for cellular structures, which eventually can 
be extended to include the pinning and entanglement 
needed for avalanches and hardness, and the slip sys- 
tems and statistically stored dislocations needed for re- 
alistic descriptions of texture evolution and cross-slip Q . 
Our model gives the elegant, continuum explanation for 
the formation and evolution of cellular dislocation struc- 
tures, ft exhibits both the experimentally observed frac- 
tal structures and scaling collapses hitherto thought in- 
compatible. Finally, it provides the fundamental distinc- 
tion between cell walls and grain boundaries; cell walls 
are intrinsically branched in a fractal fashion. 

Within a continuum theory of dislocation dynam- 
ics [13, incorporating only elastic self- interactions 




FIG. 1: Theoretical and experimental dislocation frac- 
tal morphologies (color online). Top: Simulated fractal cell 
wall pattern after uniaxial strain of e zz = 4/3o- (a) Disloca- 
tion density plot; (b) Local orientation map. Bottom: TEM 
micrographs taken from: (c) a Cu single crystal p} after [100] 
tensile deformation to a stress of 75.6 MPa and (d) an Al 
single crystal following compression to e — 0.6 Q, respec- 
tively. Gray scales have been adjusted to facilitate visual 
comparisons. Note the striking morphological similarity be- 
tween theory and experiment. 



with a minimally modified gradient dynamics, we study 
the relaxation of a smoothly deformed crystal and its sub- 
sequent evolution under external strain (Fig. QJi and b) . 
When climb is allowed, we find the distortion neatly 
evolves into a stress-free collection of grain boundaries. 
When climb is forbidden, cell wall structures evolve with 
power-law correlations and self-similarity - providing a 
clear morphological distinction between cell walls and 
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grain boundaries, a tangible model for the experimen- 
tally observed fractal structures UHl], and an alterna- 
tive to those that predict microstructure without a wide 
range of scales [12|. The resulting morphology, how- 
ever, is self-similar only after rescaling both space and 
amplitude. Performing the experimentalist's analysis of 
the misorientations and cell size distributions 0-0 yields 
good agreement with the observed scaling form (albeit 
with significantly different scaling functions and expo- 
nents). By directly exhibiting key features of the emer- 
gent experimental behavior in a continuum, deterministic 
dislocation density theory, our simulations pose a chal- 
lenge to theories based on stochasticity in the contin- 
uum laws 0, [l^ or in the splittings and rotations of the 
macroscopic cells HH. Can these stochastic theories 
describe our chaotic dynamics after coarse graining? 

Our order parameter is the plastic distortion tensor 
f3 p . Together with the resulting elastic distortion /3 E 
derivable from /3 P via the long-range fields of the dis- 
locations [Hj], f3 F both gives the deformation u of the 
material (through diUj = (3® + /??•) and gives a three- 



index variant of the Nye dislocation density tensor [16 1 
Pijk(x) = djPfk — di(3j k (defining the flux of disloca- 
tions with Burgers vector along the coordinate axis 
through the infinitesimal surface element along e; and 
ej). (3 F thus fully specifies the dislocation wall mor- 
phologies, the crystal rotation (the Rodrigues vector A 
giving the axis and angle of rotation), and the stress 
field a (the external load plus the long-range stresses 
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from the dislocations, given by a kernel 
aff +jK im {v-T')p kl {v')dv') 

Following Roy and Acharya [10j, we assume the flow 
of pijfe(x) is characterized by a single velocity v(x). Al- 
lowing both climb and glide, we can take the velocity v 
to be proportional to the Pcach-Kohlcr force J- on the 
entire population of dislocations times a mobility -D(|/o|) 
v a = D(\p\)T a = D(\p\)p ast a st , where a is the stress; 
we then define dj3 p i j/dt = Jjj = v a p a ij- (This provides 
the same equation of motion derived later by Limkumn- 
erd and Sethna [H|].) To remove dislocation climb (mass 
transport via frozen-out vacancy diffusion), we must set 
the trace of the volume change J a = 0, suggesting a 
dynamics which moves only the traceless portion of the 
dislocation density: 
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In this case, to guarantee that energy monotonically 
decreases we are led to choose the velocity based on 
the Peach-Kohler force on this traceless part v a = 
D(\p\)(past - 8 st pabb/3)crst, making the rate of chan ge o f 
the energy density the negative of a perfect square Il8|. 
(This differs from our earlier glide-only formulation jll|.) 
To ensure that the velocity is proportional to the 
force per dislocation, we choose Z?(|p|) = l/|p| = 
V y/PijkPijk/2. Our theory does not incorporate ef- 
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FIG. 2: Scaling of the correlation function (1024 2 sim- 
ulations). The trace of the orientation-orientation correlation 
function Ct 3 {R) = ((Ai(x) - A,(x + r))(Aj(x) - A^x + r))) 
is averaged over all pairs of points at distance |r| = R. No- 
tice that the simulation allowing climb has C^(R) ~ R as 
expected for non-fractal grain boundaries. Notice that the 
glide-only simulations show C^(R) ~ R 2 ~ v with n ~ 0.5, in- 
dicating a fractal, self-similar cell structure, albeit cut off by 
lattice and system size effects. 



TABLE I: Critical exponents measured for different cor- 
relation functions. GO: Glide Only; CG: Climb&Glide; ST 
Scaling Theory pi]. 



Correlation functions 



GO 



CG 



ST 



C£(r) =<£JA l (r)-A 1 (0)] 2 ) 
C*(r) =<[p^(0K-(r)]> 



1.5 ±0.1 
0.4 ±0.1 



1.1 ± 0.1 2-T) 
0.9 ±0.3 rj 



fects of disorder, dislocation pinning, entanglement, glide 
planes, crystalline anisotropy, or geometrically unneces- 
sary dislocations. It is designed to provide the simplest 
framework for understanding dislocation morphologies 
on this mesoscale. 

Our simulations show a close analogy to those of tur- 
bulent flows. As in three-dimensional turbulence, defect 
structures mediate intermittent transfer of morphology 
to short length scales. (Unlike two-dimensional turbu- 
lence, we find no evidence of an inverse cascade - our 
simulations develop structure only at scales less than or 
equal to the initial correlation length of the deformation 
field.) As conjectured [HI for the infinite-Reynolds num- 
ber Euler equations, our simulations develop singularities 



in finite time jll| . It is unclear whether our physically 
motivated equations have weak solutions; our simulations 
exhibit statistical convergence, but the solutions continue 
to depend on the lattice cutoff (or on the magnitude of 
the artificial diffusion added to remove lattice effects) in 



the continuum limit 181. Since our simulations exhibit 



structure down to the smallest scales, we conjecture that 
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FIG. 3: Self-similarity in real space (color online). Each 
frame represents the lower left-hand quarter of the previous 
frame. Frame (a) is a 1024 x 1024 simulation; (b), (c), and 
(d) are thus of length L = 512, 256, and 128. All are rescaled 
in amplitude by (L/L )~ v/2 with 77 = 0.5 (see Fig.[2]and Ta- 
ble [J). The scale is logarithmic with a range of almost 10 7 . 
Notice the statistical self-similarity. Other regions, when ex- 
panded, can show larger differences between scales, reflecting 
the macroscopic inhomogeneity of the dislocation density. 



this is a kind of sensitive dependence on initial conditions 
- but here amplified not by passage of time, but by pas- 
sage through length scales. Since the physical system is 
cut off by the atomic scale, we may proceed even though 
our equations are in some sense unrenormalizable in the 
ultraviolet. 

Wc simulate systems of spatial extent L in two dimen- 
sions with periodic boundary conditions; our deforma- 
tions, rotations, strains, and dislocations are fully three- 
dimensional. The initial plastic distortion field /3 P is a 
Gaussian random field with decay length L/5 and initial 
amplitude /3q = 1. We apply a second order central up- 
wind scheme designed for Hamilton- Jacobi equations |20| 
on a finite difference grid. The unstrained simulations 
presented are at late time, where the elastic energy den- 
sity is small and smoothly decreasing to zero, (see Sup- 
plementary Movies 1 and 2). The strained simulations 
in Fig. HI (see Supplementary Movie 3), have uniaxial 
strain in the out-of-plane direction, which is increased 
by adjusting the external stress a zz (t) to hold e(t) fixed. 
The strain rate is e = 0.05/?q. 

Figure [2] shows the orientation-orientation correlation 
function. Here we see that the cellular (climb-free) struc- 
tures have non-trivial power-law scaling, but we see non- 
fractal behavior in the grain boundary morphology allow- 



ing climb. In Supplementary Movie 2, the complex struc- 
ture of cell walls (climb-free) shows a few primary large- 
angle boundaries with high dislocation density and many 
low-angle sub-boundaries, leading to fuzzy cell walls 
that are qualitatively different from the grain boundaries 
(climb & glide, seen in Supplementary Movie 1). Table U 
includes also the correlation function of the total dislo- 
cation density; one can show [HI], if the elastic strain 
is zero (El, that C(r) = -d 2 C A (r) - d^C^r), so 
CA(r) ~ |r| a tells us that C p (r) - |r| a " 2 , implying the 
exponent relation a = 2 — r\ in the last column of Ta- 
ble |U The scaling for the correlation function for the 
total plastic distortion /3 F is not as convincing [3] ■ Both 
are consistent with a rcnormalization-group transforma- 
tion that rescales the dislocation density by a factor of 
when it rescales the length scale by a factor of b. 
Figure [3] gives a real-space renormalization-group illus- 
tration of this self-similarity; the cell walls form a self- 
similar, hierarchical structure. 

Can we reproduce the experimental fractal character- 
ization of cell boundaries? Box-counting applied to the 
dislocation density (as in Fig. QJ,) gives dimensions that 
depend strongly on the amplitude cutoff (the disloca- 
tion density is self-similar, not a simple fractal). If we 
first decompose our simulation into cells as in Fig. |4p, 
and apply box-counting to the resulting cell boundaries, 
we obtain a fractal dimension of around 1.5 over about 
a decade [laL compared to the experimental values of 
1.64 — 1.79 [3j. Such a measurement, however, ignores 
the important variation of wall misorientations with scale 
(capturing the spatial scaling but missing the amplitude 
scaling) . 

Can we reconcile our self-similar cell morphologies 
with the experimental analyses of Hughes and collabo- 
rators Using our boundary- pruning algorithm to 
identify cell walls, Fig. Hh and d show the cell size and 
misoricntation distributions extracted from an ensemble 
of initial conditions. The misoricntation distribution we 
find is clearly more scale free (power-law) than that seen 
experimentally. Under external strain, we do observe the 
experimental cell structure refinement (Fig. |4p) , and we 
find the experimental scaling collapse of the cell-size and 
misoricntation distributions (Fig. Hfc and d) and the ob- 
served power-law scaling of the mean size and angle with 
external strain (Fig.[4j: and f), albeit with different scal- 
ing functions and power-laws than those seen in experi- 
ments 

Because we ignore slip systems, spatial anisotropy, and 
immobile and geometrically unnecessary dislocations, we 
cannot pretend to reflect real materials. But by distill- 
ing these features out of the analysis, we have perhaps 
elucidated the fundamental differences between cell walls 
and grain boundaries, and provided a new example of 
non-equilibrium scale invariance. 

We would like to thank S. Limkumncrd, P. Dawson, 
M. Miller, A. Vladimirsky, R. LeVeque, E. Siggia, and S. 
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FIG. 4: Cellular structures under strain: size and misorientation distributions (color online), (a) An unstrained 
state formed by relaxing a random deformation, decomposed into cells determined by our boundary pruning method: we sys- 
tematically remove boundaries in order of their average misorientation angle, and then prune cells based on their perimeter/area 
ratio and misorientation angle (see Supplementary Movie 4). Boundaries below a threshold root-mean-square misorientation 
9 C = 0.015/?o are removed, (b) The final state after a strain of e zz = 4/3o is applied; notice the cell refinement to shorter length 
scales, (c) The cell size distribution (square root of area), scaled by the mean cell size and weighted by the area, at various 
external strains, (d) The misorientation angle distribution, weighted by cell boundary length, scaled by the mean. For each 
curve, data starts at C . This distribution appears to be closer to a power-law (inset) than the experimental distributions (solid 
curves 0-@]). (e,f) Mean cell size D av and misorientation angle Q av as functions of external strain. We find these same power- 
laws D av ~ e ~°' 26±014 and 6 av ~ e °-26±o.04^ w j^j 1 errors reflecting over a range of 8 C and for a variety of pruning algorithms and 
weighting functions. Notice that the product D av O a v is approximately constant, as observed experimentally 0]. The power-law 
dependence e 0,3 is weaker than the powers e 1 / 2 and e 2 ^ 3 observed experimentally for incidental dislocation boundaries and 
geometrically necessary boundaries, respectively. 
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